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Abstract 

We present dispersion-free wavepacket solutions to the Klein-Gordon equation, with the only free pa- 
rameter being the wavepacket velocity v . The ^'v are eigenvectors of a velocity operator v with commuting 
components, which is symmetric in a certain scalar product space. We show that v corresponds to a classical 
I generator v(x, k) which with its conjugate z(x, k) may be obtained from x, k by a canonical transformation. 

o 

■ 1 Introduction 

o 

I A particle with momentum [ra-yv) in the z direction is usually represented by the plane wave 

ip = A exp[-i ■mj{t - vz)] , |w| < 1, 7 = (1 - v'^y^^'^ 

^ solution of the Klein-Gordon (KG) equation 

^7^: (d^ -V^ + m^)il; = (1) 

^ I (we use natural units with % = c = 1 ). The plane wave is spread over all space, contrary to our usual experi- 

^ \ ence. We explore a class of solutions to the KG equation that are localised, their magnitude approximately 

inversely proportional to the distance from the centre of the wavepacket; they travel without dispersion, and 
, they are eigenf unctions of a particular velocity operator found below. The aim of this paper is to see if one 

, can replace the usual plane wave basis states with those discussed here. 

^ ■ ..... 

, The KG equation admits the following solutions (considering the wavepacket velocity v to be in the z 

. direction for the moment) 

O : vl;„(x, t) = exp[-z m^\t - vz)] , \v\ < 1, j ^ {1 - v')-'^' 

' = jo{'mj\v\r) exp[—imj'^{t — vz)], (2) 



. where r = \J + 2/^ + 7^(2 — wt)^ , 

H : 

and jo is the spherical Bessel function of order zero. The wavepacket (2) can be thought of as a plane 
wave of wavelength (27r/TO7^|u|) , bounded by the envelope function i^(m'y\v\r) which has an extent of order 
one wavelength in the z direction (for non-relativistic velocities the envelope factor jo(?n7|w|f) is nearly 
spherically symmetric about the point (0,0, ut) ). The wavepacket moves with velocity v in the z direction 
without dispersion. The ^'^ solutions are a special case of the wavepackets derived by MacKinnon. [1, 2] 
That (2) satisfies the KG equation follows trivially by applying the Lorentz coordinate transformation 

z 7(2 — ft), X ~* X , y ^ y , t ^ ^{t — vz) (3) 

to the stationary KG solution 

■01) = ia{^l\v\r) exp[— iTO7<] . (4) 

The Lorentz transformation is such that it exactly cancels out the negative z momentum components in the 
superposition making up (4), as will be seen from the Fourier transform of discussed in the next section. 
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Further solutions ^i^y of the KG equation with spin I can be generated from : 

^i,v = {dx + idyY'^y = ji{m-f\v\r)Yii{e,(j)) exp[-imj^{t - vz)] (5) 

where Yii{9, cf>) (x (x + iy/rf is the spherical harmonic. For wavepackets travelling in an arbitrary direction 
with velocity v then the '^y of (2) is 

^'v = jo("^7|v|r) exp[— im7^(f — v • x)] (6) 

where now 

f = ^Jx\-\- 7^(x|| — vi)2 , x_L = X — (v • x) V , X|| = (v • x) V . 
The (time dependant) distance f is zero at the wavcpackct centre {t, vt) . The \E'v satisfy the identity 

i (St + V • V ) *v = TO *v (7) 
noting that the operator {dt+v- V) commutes with f in (6). 

We are curious as to whether ^'v, ^v' are orthogonal, the natural approach to this question is to 
enquire if is an eigenfunction of some operator v which has commuting components. It turns out that 
this operator is easier to spot in Fourier transform space. 



2 The Fourier treinsform of ^'v 



Positive energy solutions of the KG equation are of the form (see for example Schweber [3]) 
V'(x, t) = J <^(fcO, k) e-*'=°*+*-5(fc°' - - m^) 0{k°) d'^k , 



(27r)3/2 
1 

(27r)3/2 



0(k) e- 



dk 



fc = |k| , 



<A(k) 



VrrP 



(x). 



We will call this (^(k) the momentum space counterpart to the il>{x,t) wavefunction. In the first line 
above the step function 9{k^) ensures that the energy k^ = +y/m? + is positive, and IF is the usual 
3-D Fourier transform. As ^(k) is a Lorentz scalar (due to the delta function), if we can find the (^(k) 
corresponding to the stationary solution (4) it will be straightforward to find the ^(k) corresponding to (2). 
Starting with the formula for the 3-D Fourier transform of spherically symmetric functions (j}{k) , which is 
[<j)(k)\ (x) = i / (j){k) sm{kr) k dk , then 

sin(m7|u|r) 



TT 5{k — m'^\v\) 
2 {rwyvY 



(x) 



m'y\v\r 



\ jo(TO7|w|r) 



or using the delta function identity 5{k — m7|t;|) = \v \ 5{^/rn? + fc^ — rwy) 

/7 ^(Vm^ + fc2 - m7) I 6-*^"'+*'* 



— jo(m7|u|r) exp[— im7t] 



m'y\v\ J ^rn? + fc^ 

which RHS is (4), so that the quantity in braces is the momentum space counterpart (/)(k) to (4). Next we 
substitute into the 



the transformations 



^3 7fc3 — 7f v m2 + fc2 , ki 
sj m? + fc2 — » 7-\/ TV? + fc2 + 7?;A:3 

2 



(9) 



corresponding to the Lorentz transformations (3), then we find 



TT 5{'')^/rn? + k'^ — m-y — ^vks) 
2 mjlvl 
TT S{-\/m'^ + k"^ — m — vks) 
2 



and so the (x, t) of (6) is 



*^(x,i) =:rk 0„(k) 



, — i^/m-^+k'-^t 
TT 1 



(x) 



(x) 



2 m72|t;| Vm^ + fc^ 

and the quantity in square brackets is the desired Fourier transform of "iy which we will call av(k) : 

1 



(10) 



a^(k) = 



■ 6{\/ m? + fc^ — m — i'A;3) — 



2 m72|t;| Vm^ + fc^ 

The delta function 5{^/'m? + fc^ — m — wfcs) within av(k) is non-zero on the surface 

kl + kl + 7"^[fc3 - mv^'^f = m^w^7^ , 

which is an ellipsoid tangential to the k^ = {] plane at the origin, with centre (0,0, nivj'^) and elongated in 
the /c3 direction. As already mentioned, the Q;t,(k) has no negative k^ component (the ellipsoid surface has 
fcs > 0). For the general case when the velocity v is in an arbitrary direction then av(k) is 



av(k) = 

Note that due to the delta function 



TT 1 



2 m7^|t; 



6{\/ m? + fc2 — m — V • k) 



^—i^/m'-^+k'-^t 



(11) 



-\/ + fc2 = m + V • k 

which corresponds to (7), and the time evolution operator is e~*^"*^+*'^* = e~'(™+"^''')* where the e"'"^'"** 
factor generates the translation vt as expected. 
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3 The velocity operator 

Our goal is to find the velocity operator v in Fourier transform space such that vav(k) = va(k) . Once we 
have found the operator at time t = then the operator vj for time t will be 

At time t = Q then Q!v(k) is (disregarding a factor) 

a(k) = ^(Vm2 + fc2-m -vJ^ ^^2) 

Clearly the operator v cannot be simply multiplication by a function of k as the delta function is a 2-D 
surface, not a point. Differentiating av(k) with respect to k yields a term v multiplied by the derivative 
of the delta function. We can avoid this complication by considering first the 'anti-derivative' of the delta 
function which is a unit step function: we define (x) = — 1/2, x < 0; E (x) = 1/2, x > . Then 

E{\/ m? + P — m — V • k) 

has the value —1/2 within the delta function ellipsoid, and is 1/2 outside. Then 

k 



Vfe E{y/m? + A;2 - m - v • k) = ( - v) 5{\/n^'+~k^ - m - v • k) 

= (k - y/m^+ k^v)a{k) (13) 

and 

(k • Vfe) E{^/m'^ + fc2 - m - V • k) = - (v • k) + fc^^ a(k) 



= m + A;2 — m j a(k) 

so that 

£;(Vm2 + A;2 _ m - v • k) = m (k • Vfe)"^ (Vm2 + fc2 _ ^) Q;(k) . (14) 
Now we substitute this last result into (13) obtaining 

|k - m Vfe(k- Vfe)"^ (\/m2 + F - to) | a(k) = \/to2 + fc2 v a(k) 



or 



Vfe a(k) = I ^^== - m ^=L== Vfe(k • Vfe)-^ (Vm2 + fc2 - m ) | a(k) = v a(k) . (15) 

The operator Vfe(k • Vfe)~^ within Vfe is defined as follows: 

Vfe(k-Vfe)-'= Vk{kdhr^= {-iVkk){-idkk)-'^= {-iVkk)D^'^ (16) 

where Dk is the dilation operator 

Dk = -i{dkk) 

which commutes with any operator of homogeneity degree zero (i.e. any operator invariant under a dilation 
of k). And 

D-V(k) = i{dkk)-^(P{k) = ^[J^ <^(Ak)dA - ^ 0(Ak)dA] (17) 
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which is a right and left inverse of Df- . The operator from (16) which is (— i ^kk)Df, ^ i = Z)^ ^ {—i^k) can 
be shown to be symmetric in the usual 3-D inner product space, as 

[^fc '(-*V,)]t = (-zVfe)t (Z)-i)t = (-iV,)A:i^^ii = D^\-iVu) . (18) 
From (15,18) it follows that is symmetric in the scalar product space 



ai(k),a2(k)) = (ai(k) I {^/w? + k"^ - m) \/w? + fc^ I a2(k) ) (19) 



instead of the usual Klein-Gordon scalar product space ^Q!i(k) | \/ m? + k'^ | a2(k)^ . As the components of 
Vfe which are 

{Vkl , Vk2 , Vk3 } (20) 

have a common eigenvector a(k) , it also follows that these components commute: we explicitly calculate the 
commutator [vki , Vk2\ in the Appendix and check that it is zero. It follows that eigenvectors with different 
eigenvalues of Vfe are orthogonal to each other with respect to the scalar product space (19), i.e if we label 
(k) with the property Vfeavi = victvi etc, then 

(avi (k) , (k))^ = , (vi 7^ V2) . 

In configuration space the scalar product corresponding to (19) is 

(JVi(x) \ -d^-imdt\ M^)) = \ j [^tV-t (x)] [{dt + i m)V2(x)] d^x + s.c. (21) 

where s.c. stands for the symmetric conjugate term. The operator corresponding to is v = TwkT~^ , 
found by substituting 

k— >— iV, Vfe — > — ix. 



(k- Vfe)^-(V-x) = -^(a.r)r2 =-i—Dr'' 



where D is the dilation operator D = {—idrr} into (15), obtaining 

V = —i , y = - im , — ^ 'x.\-^(iD)~^r'^] (\J m? — V^ — ml 
Vm2 - V2 Vm2 - V2 L^2V ; jvv > 

= {idt)-^(^ - iV - mx[l(£))-V] {idt-m)) 

where {D)~^ is defined as the space counterpart to (17). Recalling the ^'v of (6): 

= jo("T'7|v|f) exp[— im72(t — v • x)] 

then the result (22) implies that at time t = Q 

(V + vat)*v(x)= imx[-^(£))-V2] (iat-m)*v(x), 

which may be verified by explicit calculation, each side having the value 

— m7|t;| ^ ji(TO7|v|f) exp[— im72(t — v • x)] . 

where ji is the spherical Bessel function of order one. 



(22) 



5 



4 A canonical transformation from variables (x, k) — > (z, v) 



The fact that the components of the velocity operator commute leads to inquire whether there is a canonical 
transformation underlying this property. Classically the v operator of (22) corresponds to the generator 

1 



v(x,k) = 



k — m 



(V + fc2 — m ) 
+ fc^(k • x) 



and it may be verified that the Poisson bracket 

{va , Vb} = a, b = 1,2,3 

where the Poisson bracket (P.b.) between two variables /(x, k) , g{x, k) is defined as 

{f,g} = {df/(hz) ■ (%/ak) - (df/dk) ■ [dg/dx] . 
The generator conjugate to v is 

, , , a/to^ + k^{k ■ x) , 
z(x , k) = — — r k 



which means that 



m {\/rn?- + k'^ — m) 



{Za , Vb} = Sab , {Za , Zb} = 



as may be verified by labourious calculation. So then the transformation from the variables (x, k) to new 
variables (z, v) is canonical: 



z(x,k)= ^^^St^k 



with inverse 



v(x,k) 



X Z , V = 



m (\/rn2 + k'^ — m) 
1 



\/m?+k^ 



{\/w? + — m) 
k — m — X 



(23) 



\/rp? + A;2(k • x) 



(z.v)2 



m 



f ^ 

Vl + (Z-V)2 



k(z,v)=.^^^^S^z. 



2(z.v) 



(24) 



1 - (Z • V)2 

where z = z/2;. The value of \v\ from (23) is not limited to be less than unity, however the canonical 
transformation is meaningless if \v\ is allowed to equal unity, because the denominator (l — (z • v)^) in (24) 
can then be zero. The Hamiltonian in the new coordinates is 

■1 + (Z-V)2^ 



V + fc2 = m ^ — 



(z • v)^ 
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5 Orthogonality relation 



Wc calculate orthogonality relations between coinciding wavepackets with differing velocities in the z direc- 
tion 4*1,, vPt,' , or equivalently between their momentum space counterparts a^, a^i , recalling from (11) that 
at time t = 

'^v = i/tt - (5 (a/ m? + A;2 — m — v k^) . 

V 2 m'y'^\v\v Tn"^ + k'^ 

We have mentioned that when v ^ v' the surface delta functions inside ay, a^i only touch at the origin, 
and so whether they are orthogonal or not will depend on whether there are any factors of k multiplying 
the delta functions. We first write the above in spherical coordinates 

1 



2 m'y'^\v\^/rnF-\-W 



(5(-\/ m? + A;^ — m — vk cos9) 
^(cos^ )• 



2 m'y'^\v\vkVmF+^ vk 
Then recalling the scalar product space (• , of (19) 

(^ay , a^i^ = (cty I (\/ + fc^ — m) \/ + k^ | av> ^ 



/ 



(V + fc^ — to) \/m2 + fc^ ajja^' d{cos8k) k'^dkdcj) 

{\/'m^ + — to) Vto^ + fc^ / — 



tt"' /■ (Vto^ + A, — to) VTO^ + fc2 -/ . 

71 7 / TJT — 5 TT^ OlCOSt' 

Y'^\vv \vv J A;-^(to^+k^) 



rnrx'^'^\vv'\vv' J k^{m?+k^) vk 

X 5(cos^ ) d(cos6)k dk 

^ v'k ' 



Wm? + k'^ — to) - / \/ TO^ + fc^ — TO \/to^ + fc^ — TO. , , „ 

n :7i zz ) k dk 



m'^'y^'y'^\vv'\vv' Jo k'^y/w? + k^ v'k vk 

/-femax (^77j2 _|_ _ j^-j / w — -^77?"+"^^ — TO . \ „ 

(^;:;^)( ) * 



m'^j'^j''^\vv'\vv' Jq Vto^ + fc2 



Note that 



5(^^^^±J^)=2TO5(fc) 



but this term is annihilated by the (Vto^ + k'^ — m) factor in the integrand which is equal to k"^ /2m as 
fc ^ . Hence we can regard {v — v') as the subject of the delta function and pull out the remaining factor 
obtaining 



, (y.y> 7 



, ^ SC-^) dk 

Jo Vto2 + fc2 VV' ' 



w m'^j^Y^\vv'\vv' Jq Vto2 + fc2 
—o-^-^S(v-v') \/to2 + A:^ 



and recalling that fcmox = 2mv^'^ 



1 + 

-TO — TO 



TO^7^ 
2^2 
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6 Coordinate transformation for „ — > '^ly 

We will show that there is a coordinate transfonriation which changes the ^I/^ of (2) 
^'„(x,t) = jo{m'y\v\^y x"^ + + 7^(-s — vt)^) e^p[—im'y'^{t — vz)] 



to 



5't,'(x, i) = joim'j'lv'l^yx^ + + -/'^ {z - v' t)^ ) c^Tp[-i mj''^ {t - v' z)] . 



From inspection of the above we sec that Lorentz transformations do not preserve the form of the ^E"^ , (as 
an example the spherical solution (3) is not of the form ). The following is the required transformation: 
(t = t' 



x = mx' 



y = 



i'\vy 



yv 



f9 /? 



(25) 



7|u| ' " j\v\ 
It is readily checked that (25) implies 

j'^{t-vz) =^'^{t' -v'z') 
j^viz - vt) = i'^v'{z' - v'il) 
as necessary. Remarkably under (25) the time coordinate remains unchanged, but the space coordinate 
transformations are specific to both the original and transformed velocities v and v' . So that if the transfor- 
mation (25) is applied to, say, ^i^i with v" ^ v , the resulting function does not in general have the form 

^ 11 • 



7 Outlook 

The wavcpackct solutions ^'v(x, t) of the KG equation arc eigenvectors of the velocity operator (22), whose 
components commute and are symmetric in the inner product space (21). These wavefunctions are interesting 
alternative basis state to the usual plane waves. 
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Appendix 



Recalling 



— m ■ 



with Dk = —i {dkk) and 



m , 



</)(Ak) dA 



we will need the following identities: 

W ' fcVi] = 



[£>^^ , A;iV2 - feVi] = 

A: n . m^fc 



k 



(Ai) 

(A2) 
(^3) 
(A4) 

(A5) 



To calculate [■ui , ■U2] we first calculate those parts of the commutator with the ^^^^ terms, which are 



m ■ 



V m2 + A;2 ' V rn? + fc2 



1 



\/m2 + fc2 ' \/m2 + fc2 



fc2 



TO 



= Z ■ 



fc2 



(fci V2 - fc2Vi) kDl^\ (Vm2 + fc2 - m ) - 



TO , „ , s Vto2 + k"^ —m 



1 



Vto2 + fc2 LVto2+ A;2 ' 



= im{kiV2 - fc2Vi) ^ 
Secondly we calculate how the second terms of commute with each other, yielding 



(a/ to2 + A;2 — to ) ^ 



i?fe^)-(v/m2+ /e2-TO)). 



v'to2 + A:2 



(A6) 



TO 



=TO 



Vto2 + fc2 -^fe 
Vto2 + fc2 



' ((-iVi) 



•\/to2 + A;2 



V to2 + fc2 



\/to2 + fc2 " ^ ' ' ' Vto2 + fc2 

(-zV2fc) + ^ , „ . (-«Vi/c) 1- 



^fe ' (-«Vi) ) ( + fc2 - to) 



= - TO^ (fciV2 - k-iSx) 



s/m? + k^ 
1 



Vto2 + fc2 Vto2 + fc23 

Adding (A6) and (A7) yields zero. 



s/w? -\- A;2 

^fc'^(Vm2+ fc2-TO) 



z2 + A;2 — to) 



(A7) 
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